Curvature effect on the interaction between folded graphitic surface and silver clusters 
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Evidence of curvature effects on the interaction and binding of silver clusters on folded graphitic 
surfaces has been shown from both experiment and theory. Density Functional Theory (DFT) 
calculations within the local density and generalized gradient approximations have been performed 
for the structural relaxation of both Ag and Ag2 on curved surfaces, showing a cross-over from 
quantum to classical behaviour. Using Lennard- Jones potential to model the interaction between a 
single cluster and the graphene surface, evidence is found for the curvature effect on the binding of 
silver nano-particles to folding graphitic surfaces. The theoretical results are compared to SEM and 
AFM images of samples obtained from pre-formed silver cluster deposition on carboneous substrates 
exhibiting anisotropic pleat structures. 
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Introduction. In the quest for nano-scale science 
and technology, interaction and diffusion of adatoms, 
molecules and nanoparticles on surfaces attract enormous 
interest due to their relevance in the construction and 
the stability of new nanoarchitectures. Most studies of 
adatom diffusion have focused on flat surfaces and it is 
only since a decade ago that the study of diffusion on 
deformed surfaces has opened up new perspectives in 
anisotropic diffusion. It has been shown that the interac- 
tion of silver atoms with carbon nanotubes is curvature 
specific [H- Both experiments and simulations pointed 
out that convex surfaces, such as those on the outside 
of carbon nanotubes, enhance binding of silver atoms. 
In contrast, our recent experiment on silver cluster de- 
position on folded graphite shows evidence of repulsive 
barriers for convex bends Q ■ 

On graphite surfaces, metal clusters of up to a few thou- 
sand atoms are known to be highly mobile. They ag- 
gregate into fractal islands @, 0] through an isotropic dif- 
fusion of clusters on flat terraces, which has been un- 
derstood using a Diffusion Limited Aggregation (DLA) 
modeljfl]. The islands with anisotropic morphologies ob- 
tained on curved surfaces reveal anisotropy in cluster dif- 
fusion, which was attributed to surface curvature. In 
fact, the weak interaction between two adjacent graphene 
sheets show evidence of pleats, observed via AFM imag- 
ing techniques. These pleats destroy the homogeneity of 
the graphite surface, which should affect cluster mobil- 
ity and thus change the island patterns @. Understand- 
ing the role of the surface curvature on cluster mobility 
becomes of fundamental interest for building new archi- 
tectures. However, nothing is known on the mechanism 
underlying the diffusion, and in particular the binding in- 
teraction between the cluster and folded graphite has not 
been studied. Two main questions need to be answered 



to fully understand and predict the final growth pattern 
from nanoparticle diffusion and aggregation on surfaces. 
Since atoms and clusters diffuse with opposite behaviours 
on convex curvature, at which cluster size does the tran- 
sition from atom to cluster occur? And, for clusters, how 
does the cluster-surface interaction depend on the sur- 
face curvature? 

In this letter, we show evidence of a curvature effect on 
the interaction between silver nano-particles and folding 
graphitic surfaces, both experimentally and theoretically. 
We show that the islands grown from silver nanocluster 
deposition on graphite remain fractal on long range un- 
dulated surfaces, as it is for fiat terraces, whereas those 
grown on more concave region of the pleats are elongated 
structures that are located in the valley of the pleat. Sev- 
eral levels of theoretical methods, from quantum mechan- 
ics to continuum model, have been used to address var- 
ious aspects of the problem. Using Density Functional 
Theory we note that the binding energy of an atom and 
a dimer on a perfect carbon nanotube behave in oppo- 
site ways. For larger clusters, using a Lennard-Jones po- 
tential to model the interaction between a single cluster 
and the graphitic surface, we provide quantitative infor- 
mation on binding energy changes as function of surface 
curvature. The theoretical results fully explain the exper- 
imental observations on cluster deposition on carboneus 
substrates. 

Experiment. In our experiments, a distribution of neu- 
tral silver clusters with a mean diameter of 3 nm and half- 
width at half maximum of 0.5 nm is deposited at thermal 
energy on folded graphite. The low impact energy of 0.05 
eV/atom, as compared to the Ag-Ag 0] binding energy 
1.2 eV, makes the fragmentation of the impinging clus- 
ters unlikely. They diffuse on the surface as a whole and 
grow into islands. The island morphology of the sam- 



pies is separately analyzed by SEM and tapping AFM. 
On graphite terraces the clusters aggregate to fractal is- 
lands anchored on point defects as described before 0, Q] . 
In this paper we focus on cluster deposition on curved 
graphitic surfaces. The pleats are typically many \im 
long and resemble macroscopic drapery even though 
the graphite is mostly facetted like graphite polyhedral 
crystals as drawn in Fig. [Ha). We focus on three- faced 
pleats, whose width ranges from 100 nm to a few microns. 
Atomic force microscopy (AFM) reveals that their height 
ranges from 10 to 35 nm. and the bend angles 9 of the 
pleats is quite low - typically range between 3 and 35 
degrees. Although the angles of the pleats are quite well- 
defined from their AFM profiles, the estimated values of 
the radius of curvature are less accurate and are limited 
by the size of the AFM tip - typically — 10 nm. Figure 
QJb) and QJc) show typical SEM and AFM images of a 
graphite pleat of 250 nm width after silver cluster depo- 
sition. As described earlier the clusters aggregate to 
linear islands in the concave bends of the pleat, in which 
they are trapped, but not demobilized as on point defects 
or step edges. We reported further, that convexly bent 
areas act as effective repulsive barriers for the cluster dif- 
fusion, which the clusters cannot pass by their thermal 
diffusion. Such potential barriers are therefore at least 
25 meV high (thermal energy). By contrast, the islands 
grown on slightly curved graphite show isotropic fractal 
morphology that settle with the shape of the graphitic 
surface (Fig. [ljcl)). This indicates that cluster mobil- 
ity is not affected by graphite curvature of small angle 9 
equal to 3 degrees. 

DFT calculation. To compare the behaviour of atoms 
and clusters on curved surfaces, we start by calculating 
binding energy of a silver atom and dimer on curved 
graphene. We used perfect carbon nanotubes (CNTs) 
to model the curved graphene. The electronic structure 
and structural relaxation calculations were performed 
using Density Functional Theory (DFT)Q] within the 
local density and generalized gradient approximations 
as implemented in the Quantum-ESPRESSO0 pack- 
age. In our calculations, plane wave basis sets, PBE 
exchange-correlation functional[l(3] and LDA potential, 
and RRKJ ultrasoft pseudopotentials ll( have been 
employed. We used nonlinear core corrections for the Ag 
atoms, with the 3d state included in the valence. The 
use of ultrasoft pseudopotentials enabled us to use an 
energy cut-off of 32 Ry for the plane wave basis, while 
the density cut-off was taken to be 400 Ry. The Brillouin 
zone was sampled by 1 x 8 x 1 special k-points using the 
Monkhorst-Pack scheme [ijj], and a Gaussian smearing 
of width a = 0.001 Ry was used for electron occupations. 

We started with a single Ag atom, comparing the bind- 
ing energy as calculated by DFT-GGA on the concave 
and convex surface. On a (5,5) nanotube, the silver atom 
is more strongly bound to the convex surface (outside) by 




Figure 1: (Color online) (a) Schematic drawing of a pleat with 
the top width 1 and bottom width L, as well as its concave 
radius of curvature R and angle 0. (b) and (c) SEM and 
AFM images of a graphite pleat of 250 nm width after cluster 
deposition. Notice that islands grow only on the top faces and 
convex valley. 9 equals 18 and 7 degrees for the concave and 
convex bends, respectively (note that the pleat is asymmetric 
and the convex bends were not measured) (d) AFM image of 
a larger pleat, few fim width and comparable height to (c). 
In this case the angle of curvature is 3 degrees and fractals 
patterns are not affected by the curvature. 
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Figure 2: (Color online) (a) Adsorption of Ag2 on the outer 
and inner surfaces of a (10,10) CNT, d.Ag-C is the average of 4 
nearest Ag-C distances, (b) distribution of the first 33 Ag-C 
distances, the filled box corresponds to Ag2 outside the CNT 
(left in panel a), and shaded outside the CNT (right in panel 



0.20 eV. On a (10,10) CNT, which has lower curvature, 
the two binding energies are nearly equal. We then com- 
puted Ag 2 on a (10,10) CNT (see Fig. [2{a)) and found 
that the dimer prefers inside site by ~0.3 eV 

In a recent work by Jalkanen et al., it was noted that 
the van der Waals interaction is a large component of 
the binding between silver and graphene, and GGA may 
not be an adequate approximation for the small curva- 
ture limit (l3l| . To examine this issue, we repeated the 
above calculations using LDA. Our results show that the 
Ag atom energetically prefers the convex surface of the 
(10,10) CNT by ~0.1 eV. As in the GGA calculations, 
the dimer Ag2 is more strongly bound to the concave sur- 
face of a (10,10) CNT by -0.2 eV, in contrast to the sin- 
gle atom. The unbound atom prefers the convex surface 
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where the stretched C-C bonds allow for some chemical 
binding. The dimcr, on the other hand, displays classical 
behavior, i.e. it can be modeled by a pair-wise additive 
potential. 

The concave surface provides a closer area of contact, 
as illustrated in Figure EJb) , and thus higher binding. 
This means that Ag„ clusters can be modeled classically 
for n>2. 

For single Ag atoms, quantum-mechanical effects need 
to be taken into account explicitly. A further argument 
for the classical picture for n>2 comes from the mag- 
nitude of the binding energy. To date, there is no di- 
rect experimental data on the desorption energy of a Ag 
atom or cluster on a graphite surface. Theoretical cal- 
culations give values of binding energy and bond length 
with a large error bars (a few tenth of eV). We have 
performed LDA calculations for a Ag2o cluster on a flat 
graphenc sheet. The binding energy between Ag2o and 
the graphene sheet was 1.20 eV. The clusters are rela- 
tively weakly bound, again confirming the validity of a 
classical model. 

Model calculation. We use a Lennard-Jones potential 
to model the interaction between a single cluster and the 
graphenc surface. The cluster is far enough away such 
that the fine detail of both the cluster and the graphene 
sheet can be approximated by a solid sphere and continu- 
ous sheet, respectively. We can then write the interaction 
using a standard Lennard-Jones potential [3] , scaled to 
make it dimcnsionlcss, 
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where e is the interaction energy, a is the interaction 
distance parameter, R denotes the cluster position and 
f is integrated over the infinite graphene sheet. To more 
properly model a cluster, we integrated the above ex- 
pression over the spherical volume as well. However, the 
resulting potential obtained is qualitatively no different 
from a Lennard-Jones potential with different parame- 
ters, so that step is neglected and we shall consider the 
cluster as a single entity at the sphere's origin. All dis- 
tances can simply be scaled by a, so the absolute value 
does not give quantitatively different results. We used 
this model to study the interaction between the cluster 
and a pleat in the graphene sheet. The pleat was mod- 
eled by two straight sheets at the pleat angle 9, inscribed 
by a cylinder segment of varying curvature k. to soften 
the sharp angle (see Fig. [3]). 

Figure [3] shows the potential for a 15° bend, inscribed 
with a cylinder of k = l/2 1 / 6 a. The concave bend shows 
a minimum relative to the flat pieces at the edge of the 
plot, which corresponds to stronger binding. The mini- 
mum occurs due to increased surface area for interaction 
as clusters approach the pleat. For the convex bend we 
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Figure 3: (Color online) Top: Diagram of model bend. Mid- 
dle: The potential surface, scaled by 4e for a, 6 = 15° concave 
bend. For clarity only the values around the minimum are 
colored. The lines near the bottom show the cylindrical seg- 
ment (ft = l/2cr for these plots) (red) and the straight planes 
(black) used to model the bent graphene. All distances are 
scaled by a. Bottom: Same as middle for a convex bend. 



observe the opposite, a decrease in binding due to de- 
creased nearby surface area. To make this more quanti- 
tative, we considered the binding energy relative to the 
flat surface as a function of curvature. We define the rel- 
ative binding change A = E curve /E f i at - 1, where E flat 
is the energy minimum found for a flat surface. Figure [4] 
shows A as a function of the cylinder curvature n. 

For very low curvature (large cylinder radius), the 
binding energy of both concave and convex bends ap- 
proach the flat result. In the limit of zero curvature, they 
are identical. In the limit of large curvature (small ra- 
dius), there is strong enhancement of binding for 6 > 0, 
and suppression for < 0. The inset shows that, as 
a function of 9, there is a binding energy change for a 
wide range of curvatures. Furthermore, it is approxi- 
mately symmetric about the origin - the enhancement 
corresponding to a concave bend is roughly equal to 
the suppression corresponding to a convex one. For the 
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Figure 4: (Color online) Semilog plot of A (see text) as a func- 
tion of cylinder curvature. The pleat angle 9 increases mono- 
tonically from —15° to 15°, top to bottom. Inset: Constant- 
curvature cuts as a function of 9. 



ranges of angles observed in experiment, there enhance- 
ment/suppression ranges from 2% to over 10%. 
Discussion. The calculations show that typical bends, as 
found in experiments, can either enhance or suppress the 
total interaction by up to 15% compared to a flat surface. 
Our DFT calculations indicate a total binding of around 
1 eV for Ag2o on the planar graphitic surfaces, which is 
comparable to previous theoretical studies [l5l. fli|. 

We use this value as a rough lower limit for bigger 
clusters like Agsoo, as used in our experiment. Based 
on this, the calculations show that there is a minimum 
bend angle required to trap the clusters in the bend. 
Energetically, the trapping starts between 2% and 4% 
enhanced binding (at room temperature). This corre- 
sponds to a bend angle between 3° and 5°, as long the 
radius of curvature of the bend does not exceed 10 times 
the cluster radius, or about 15 nm for Agsoo. Similarly, 
convex bends with corresponding parameters cause effec- 
tive barriers that cannot be traversed by cluster diffusion. 

The results of our model calculation explain the exper- 
imental observation shown in Fig. [TJ When clusters mov- 
ing on the graphene surface encounter a concave bend of 
a pleat, they are trapped by the enhanced binding (po- 
tential well relative to the flat surface). They now have 
to diffuse in the potential well along the bend and aggre- 
gate to linear islands. Clusters landing on top of a pleat 
that has a flat surface are unable to overcome the barrier 
at the convex bends, and are trapped on top of the pleat 
where they have to aggregate. 

By contrast, clusters moving on graphitic surfaces with 



an overall small angle of curvature do not " see" the cur- 
vature and aggregate as shown on Fig. [TJd). Our DFT 
calculations show that the bonding of single silver atoms 
is oppositely influenced by the surface curvature. Con- 
sequently, the manner in which bends in graphite affect 
atoms and clusters is qualitatively different. The DFT 
calculations further show that silver dimers already be- 
have classically and cluster-like. As the cluster size is in- 
creased further, only the total binding, and therefore the 
effective wells or barriers, caused by bends in graphitic 
pleats increases. For clusters with 500 silver atoms the 
potential wells and barriers are high enough to be used 
as effective guides or traps for cluster diffusion and ag- 
gregation, opening new routes for controlling anisotropic 
diffusion at nanometer scale. 
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